EFI-94-42 
hep-th/9409128 



Wilson Loops, Winding modes and Domain Walls 
in Finite Temperature QCD 



Joshua Boorstein and David Kutasov 

Enrico Fermi Institute 
and Department of Physics 
University of Chicago 
Chicago, IL 60637, USA 

We discuss the effective action for Polyakov- Wilson loops winding around compact Eu- 
clidean time, which serve as order parameters for the finite temperature deconfinement 
transition in SU{N) Yang-Mills gauge theory. We then apply our results to the study of 
the high temperature continuation of the confining phase, and to the analysis of certain 
Zn domain walls that have been argued to play a role in cosmology. We argue that the 
free energy of these walls is much larger than previously thought. 
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1. Introduction 

There are numerous reasons to study SU (N) gauge theory at finite temperature. In 
particular, detailed understanding of the SU (3) theory should be useful in describing the 
physics of QCD in the early universe, and properties of the quark-gluon plasma [|1| (the 
high temperature phase) are important for relativistic heavy ion collisions. More formally, 
one is interested in understanding the deconfinement transition at large and study its 
relation to the ideas of Hagedorn and to string theory . This is especially interesting 
since the perturbative gauge field description, suitable at high temperature and the string 
one, valid at low temperature (in the confining phase), should be complimentary to each 
other. 

As we will briefiy review later, the free energy e~^^ = Tre~^^ and other physical 
quantities are given at finite temperature by a path integral over gauge fields living on 
the Euclidean manifold x with Euclidean time identified with + /3; gauge 
fields are periodic: A^(x° + /3,x) = A^(x°,x), while quarks (if present) are anti-periodic. 
There is a nice correspondence between the phase structure of SU{N) gauge fields and Zn 
spin systems 0]. The high temperature (deconfined) phase in gauge theory corresponds 
to the ordered (low temperature) phase of the spin system and vice versa. The scalar 
gauge invariant order parameters which capture the dynamics of gauge fields are time-like 
Polyakov - Wilson loops: 

(with the trace in the fundamental representation of SU{N)). Only — 1 of the Wn 
are independent. Their correlation functions are of utmost importance; e.g. the quark- 
antiquark free energy Fq^q defined by: 

g-/3F,,,(x,y) _ (^Wi{^)W-i{y)) (1.2) 

measures the free energy of a system with a static quark at x and an antiquark at y. Higher 
Wn are related to higher representations of SU{N). In the confining phase {Wn) = and 
g(x) ~ |x| as |x| oo, while in the deconfined phase {Wn) 7^ and ^(x) ~ constant 
as Ixl 00. 
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To focus on the dynamics of the Wn one may integrate out the other degrees of freedom 
and study an effective action of the general form: 

5eff =iV2 / c/Wy5^M^„(x)W^_„(y)G'(f)(x-y) 

- n 

+ / d^Xid^XaC^^Xa J] VT^, (xi)Vr^, (x2)W^-ni-n2 (Xa)^^^,^,^, (Xi, X2, X3) + ■ ■ ■ . 

ni,n2 -' 

(1.3) 

The kernels G^'^\G^^\.. . summarize the dynamics, and are the object of this paper. Eq. 
( |T73|) is the analog of a Landau- Ginzburg description of the spin system, and can be used 
in the usual way to study the behavior of Fq^q (|1 . 2|) and other observables. At large the 
structure is similar to string theory. Terms of high order in Wn are suppressed by powers 
of (/string = In particular, it is very interesting to compare the inverse propagator for 

Wilson loops G^"^^ to the corresponding quantity in conventional string theories; knowledge 
of G^"^^ is sufficient to deduce the value of the large deconfinement (Hagedorn) transition 
{(^h) and other physical properties, like the string tension as a function of temperature. 
The immediate motivation for our work is related to three recent ideas: 

1) In 0, J. Polchinski proposed to study properties of the confining phase (essentially 
properties of G^"^^ ( p..3|) ) of four dimensional Yang-Mills theory using perturbation theory, 
valid at high temperature, by "analytically continuing" the confining phase to that regime. 
These arguments were since generalized to other cases 0; there have also been attempts 
to match the resulting behavior of G^'^^ to particular string theories (strings with Dirichlet 
boundaries [H, and rigid strings ||]). As pointed out by Polchinski it is surprising 
to obtain properties of an essentially non-perturbative object (the confining phase) from 
perturbation theory. Our study of (|1.3| ) will allow us to clarify this issue somewhat. 

2) It has been proposed that at high temperature, when one knows that the minima 



of S'eff are at (VFn) = e ^ \ j=0,l,...,A^ — 1, the system may possess domain walls sep- 
arating regions in space with different (VFn); these domain walls may be of cosmological 
interest. The interface tension (energy per unit area), a, of such domain walls has been 
argued to be perturbatively calculable and to go like a ~ -^^^ (with T the temperature 
and g the running gauge coupling) . If valid, these arguments could suggest interesting new 



phenomena that may have occurred in the early universe [jT2[. In addition, these domain 
walls would correspond to non - perturbative effects of order exp(— l/gr) in (finite temper- 
ature) gauge theory, which would be of more general theoretical interest [jl^. However, 



various objections to this scenario have been raised [14], [15]. In particular, it has been 
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claimed that infrared (IR) divergences may lead to subtleties in the arguments of |11 
although the leading behavior a ~ was argued to be safe. We will reconsider some of 
these issues below. 



3) In 0, |T^ it has been argued that two dimensional QCD coupled to adjoint "quarks" may 
be a good toy model of the relation between large gauge theory and strings, exhibiting 
a non - trivial spectrum of "Regge trajectories" and a large deconfinement transition. 
It is of interest to develop techniques that would allow calculations of properties of this 



Hagedorn transition. The results of this investigation will be reported separately |T8 

In this paper we are going to discuss the structure of the effective action S'eff, presenting 
techniques to evaluate G^'^^ perturbatively in the gauge coupling. The plan of the paper 
is as follows. Section 2 is devoted to a brief review of finite temperature gauge theory, 
presented mainly to set the notations and specify the necessary calculations. In section 3 
we describe the calculation of G^"^^ to one-loop order and outline the calculations of higher 
order kernels {G^^\ ■ ■ ■). The perturbative analysis is seen to be infrared finite in a certain 
region. We discuss separately the contributions to G^'^'^ of gluons and other kinds of adjoint 
matter. 

In sections 4, 5 we discuss the lessons learned from the one-loop calculations of section 
3 for the two problems mentioned above, of the high temperature limit of the confining 
phase, and the interface tension of domain walls. We show that infrared divergences do 
not alter the analysis of 0, (contrary to recent claims), but unfortunately one does 
not appear to be able to learn much about properties of the confining phase from this 
analysis. On the other hand, for the domain wall problem infrared issues are found to play 
a crucial role, and in fact alter the qualitative behavior of the domain wall energy per unit 
area a [jTl|, from a ~ ^ to a ~ p-. Section 6 contains a summary of our conclusions and 



necessary future work. 



2. General Formalism 



It is customary 0, to calculate quantities like 

= Tre"^-^ (2.1) 

in the Aq = gauge. Ignoring matter fields for simplicity, we have for the QCD Lagrangian 
C = ^TrF^ the Hamiltonian: 

n=lj rfW(E")' + ^(B«)2] (2.2) 



where E",B" are the color electric and magnetic fields respectively; E",A" are canoni- 
cally conjugate. The physical Hilbert space is spanned by \A{x)) satisfying Gauss' Law 
constraint 

D • E|phys) = 0. (2.3) 
This constraint may be enforced via a Lagrange multiplier: 

Z = J VA{x){A\e-f^^P\A) (2.4) 

where P = J T>T{x)e^ -f xTTDr{x)-'E{x) p],Qjg(,^g solutions of Gauss' law (|2.3|). Standard 
Feynman path integral manipulations then lead to the expression: 

Z = J [VA^{x)]e~^ -fo d^'Jd'^^<. (2.5) 

where F has been renamed Aq and the gauge fields have the periodicity properties: 

A^(x° + /3,x) = A^(xO,x). (2.6) 

One can also study more sophisticated questions like what is the free energy in the 
presence of sources. For that one needs to generalize the Gauss' law constraint (|2.3|) to 



D • E = p. Repeating the previous discussion one finds Q , |jT^ that the free energy with 
static quarks at positions xi,...,Xn and antiquarks at yi,...,yn is given in terms of 
Euclidean time-like Wilson loops (|1.1|) by: 

n 

g-/3F(x„...,x„,y„...,y„) ^ ( JJ VTi (Xi)T^-l (yi)) (2.7) 

i=l 



with the average performed in the measure ( p.5|) . We see that the dynamics of the Wn 
encoded in ( |1.3| ) summarizes many important physical properties of the theory. 

The theory (2^) is invariant under gauge transformations — > A^U -\- U~^d^U 



such that U{x'^ + /3,x) = 2t/(x°,x). The prefactor z is constrained by the periodicity of 
A^ ( |2.6| ) to lie in the center of the gauge group. For SU{N) the center is simply and 
z = , n = 0, 1,...,A^ — 1. Local gauge invariant observables are invariant under these 
aperiodic transformations, while Wn (|1 . 1| ) transforms as: 



Wn ^ z^Wn. {2.i 
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The effective action (|1.3|) must exhibit the global Z^v symmetry (|2.8|) , but this symmetry 
can be broken spontaneously. Indeed, in the confining ("disordered") phase where this 
("center") Z^r symmetry is unbroken, {Wn) = 0, while at high temperature (Wn) 7^ and 
Zjsf is spontaneously broken just as in the ordered phase of the spin models 0]. 

To calculate quantities like ( p.7| ) in the continuum, one may proceed as follows. First 
fix the gauge; ( p.5| ) , (|2.6| ) does not allow going to = gauge. A convenient gauge choice 
is 

A-,\x^,^) = jeai^)S-'. (2.9) 

Thus Aq is diagonal and independent of xq and = (mod 1). Then integrate out 

the spatial gauge fields Ai in ( ^75|) , and whatever (adjoint) matter fields are present, and 
find an effective action of the general form (|TT^) written in terms of ^a(x) or ( |1 . 1| ) : 



1 ^ 

Vr„(x) = _ ^e^^^-^W. (2.10) 



a=l 



The description in terms of ^a(x) is natural in the deconfined phase (at high tem- 
peratures) where, as we will see, the effective action is minimized for certain fixed 9a, so 
that (Wn) 7^ 0. Below the deconfinement transition the 9a are randomly distributed and 
(Wn) = 0, thus it is preferable to describe the dynamics in terms of Wn- We will use 
perturbation theory in the gauge coupling g, which in principle should be reliable at high 
temperatures (up to possible IR divergences) due to the running of the gauge coupling; 
nevertheless, we will mostly write the effective action in terms of ( |2.10| ), in the spirit 
of H ; the Wilson loops seem to be the suitable variables for discussing the Hagedorn tran- 



sition |T^, and may teach us something about high temperature string theory. To study 
the deconfinement transition and/or the confining phase one must use non-perturbative 
tools |TH]. 



3. The One-Loop Effective Action 

We will be mostly interested in two general classes of gauge systems: 
1) Two dimensional Yang- Mills coupled to adjoint (bosonic or fermionic) matter 0, |17]. 
The Lagrangian for scalar matter is 

/: = (D^cP)^ + m'<P^ + ^F^ (3.1) 
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where 

D^(j) = d^(/} + i[A^,4>]; F = doAi - diAo + i[Ao, Ai]. (3.2) 

We choose the gauge (|2.9|) ; the Faddeev-Popov determinant is det(9o +iAo). Adding to 
this the effect of integrating out Ai, ^ to one-loop order we find 

g-si-^°°p(Ao) _ [det[(ao + iAof]]'^ (det(-L>2 + m^))"* det(ao + iAo). (3.3) 

with the three determinants on the r.h.s. arising from integration over Ai, (p and the 

ghosts, respectively. The ghost contribution exactly cancels that of Ai and we are left 
with 

5,V°°P(Ao) = ^logdet(-^2 +m^). (3.4) 



For constant Aq one can use the results of [0, to evaluate SesiAo). We will discuss 
the determinant for an x dependent Aq ( |2.9| ). 

For fermionic (Majorana) adjoint matter one finds 

= '1 logdet{-DH + F^^J^^ + mH) (3.5) 

where J^'^ = |[7'^,7'^] are the Lorentz generators and = — lA^. 
2) D-dimensional pure Yang-Mills (with D>2). Thus we start with 

C = ^Fl^. (3.6) 

gZ 

Here, it is convenient to use a background field gauge and Feynman gauge for the quantum 
fields. The one-loop effective action in this gauge is 

Sl^'^^^{Ao) = \ logdet(-D2l + F^^jn - logdet(-D2) (3.7) 

with the Lorentz generators (J^j^)'"^ = ^(?7^?7^ The second term on the right hand 

side of (|3.7|) is due to the ghosts, whose contribution is minus that of a (complex) massless 
scalar in the adjoint representation. 

The one-loop effective actions ( |3.4| ), ( |3.5| ), ( p.7| ) should be evaluated for = Aq^^o 
(PI) and added to the classical action Sfi"" = ^ J rf^"^x(VAo)2. 

Most of the analysis will be done for the case of scalar adjoint matter, which (at one- 
loop) has most of the qualitative features of the other two cases. The properties special to 
the other cases will be mentioned later 0. 



^ The following two subsections are rather technical. Readers interested only in the results 
may proceed directly to eq. (3.34). 
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3.1. The one-loop effective action for scalar adjoint matter - general considerations 
It is useful to use a first quantized representation to calculate the determinant 



(see e.g ||21[] for a recent review). The effective action is given by a path integral over 
worldline trajectories x^{t),0 < t < T that wind n times around the compact xq direction: 



x(t + T)=x(f); xo{t + T) =XQ{t)+n(3. 



(3.8) 



One has also to sum over the length of the world line T G R"*" and winding number n G Z. 
Thus the one-loop effective action ( ^.4] ) can be written as: 

n=-oo J \_ Jo 

(3.9) 

with the trace taken in the adjoint representation and A/" a normalization factor. This 
form of expressing the determinant ( |3.4| ) is related to the standard Feynman diagram 
representation by Poisson resummation. It is convenient to exhibit the n dependence 
explicitly by taking xo(t) — ^ xo{t) + such that xo(i) is now periodic and the effective 
action becomes: 



n — — •-'0 ^ 



n= — oo 
T 



exp 



dt 



2 02 



1-2 , n'f3 
4^ + 4T2 



~iAo[x{t)] xoit) + 



m T 



(3.10) 



The path integral over x{f) for a general Aq ( |2.9| ) defines a complicated quantum 
mechanics problem (which is closely related to the amplitude for pair production in a 
general external electric field), but fortunately we do not need the full solution, at least 
when trying to set up the expansion (|1.3| ) in powers of Wn- As a warm up exercise, consider 
the simple case of constant Aq (or 6a (PT^)). Due to the periodicity of xq we have: 

n(3AQ, 



n£ 



dtAo (^xo{t) + 

and (replacing SU{N) by U{N), which is unimportant at large N) 

11 ^ dT 



/'l — loop 



^ (47r)2 a,b=ln=-oo 



(3.11) 
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where Cn = (^ 27rn/3 j Kd{nm[3) and K^i is a modified Bessel function. For m = 0, (i = 4 
this gives the result of (for constant Aq there is no difference between the scalar ( |3.4| ) 
and the gauge ( p.7| ) determinants), while for (i = 2 it agrees with 0. We also notice 
that the winding of the adjoint scalar around compact time n in (|3.1C1| ) is identical to 
the winding of the Wilson line n in (|1.3D (using ( |2.1C1| )). This is a general feature of the 
one-loop determinants (as we will see below) which will make things easier later. Note also 
that for constant Aq (or Wt^( |1.1[ )) the effective action ( |1.3p is exactly quadratic to one- loop 
order, i.e. G^^) = G^^^ = . . . = 0. We wiU see soon that this is not the case for non-zero 
momentum and for higher orders in the gauge coupling. 
For arbitrary Aq[x) one may proceed as follows. Write 

x{t)=x,i + Xq{t) (3.12) 

where Xd is constant, and Xq is the fluctuating quantum mechanical variable in ( |3.10| ), 
which is Gaussian with: 

«(ti)<(t2)) = -G{ti - t2)g^''; G{t + T)= G{t) = \t\ - t^/T. (3.13) 

Now substituting ( p.l2| ) in ( p.lO|) and expanding (for simplicity we take Aq to depend on 
only one of the d—1 spatial directions; no generality is lost because of rotational invariance) 
we get: 

1 1 

Aq{xci + Xq) = Ao{xci) + AQ{xci)xq + -Aq{xci)xI H h — (Xcz)Xg H . (3.14) 

^ Til . 

We will write this compactly as Ao[a;(t)] = Ao[x] + Ao[x(t)] where we use x to denote Xd 
and x{t) to denote Xd + Xq. 

Now expand (|3.10| ) in powers of Aq (or 9a) and use ( |3.13| ) and ( |3.14| ) to average over 
quantum fluctuations of the trajectory. We find that S'^g [Aq] is naturally expressed in 
terms of bilinears in quantities of the form {oj^^^ denotes the n'th derivative of 9a)'- 

1 ^ 

Vn:l,n, = -^^^)...e'=)e2™^''(--); > 1, (3.15) 

a=l 

i.e. 

cl-^^^^[v] = iv^ $:«{-}w.}<\^Sr (3-16) 
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This structure follows immediately from properties of the trace in the adjoint repre- 
sentation (for U{N)). A Wilson loop in the adjoint representation of U{N) can be written 
as a product of a Wilson loop in the fundamental and one in the anti-fundamental N 
representation, i.e. 



^1 — loop 



1 1 

2 



(47r)2 

1 1 

2 (47r)^ 
TrN exp 



E 
E 



dT 



.g 4T^-m T 



Tf + 1 



Tradj exp 



dtAo ■ \ xo + — 







g 4T ^ 



7 4 I • 

dtAo ■ I xo + — 



TrN exp 



1 A t ■ 

dtAo ■ xo + — 



(3.17) 

Now for the gauge choice ( |2.9| ), each Wilson loop in the fundamental representation can 
be written in the form (|3.15| ) by first expanding the gauge field as in (|T|) and then 
expanding the exponent in powers of 9, i.e. 



Tr^r exp 



dtAo ■ {xo{t) + Y- 



N 



N 

E 

a=l 



,27rm6'o(x) 



1 + 



2ni 
TJo 



dt6aix{t)) xoit) + 



a=l 



4n 



2 rT 



dU 



dt29a{x{ti))9a{x{t2)) [ Xo{ti) + 



,2-Kin9a{x) 



N 

E 

a=l 

47r2 rT 



1 + 



2ni 

TJo 



nl3 



Xo{t2) + ^ I + 



dt 9'^ix)x,it) Xoit) + ^] + 







dti / dt2 9'^{x)9^{x)xq{ti)xq{t2) xo{ti) -f — xo{t2) + — + 



nl3 



nl3 



+ 



W(-^ + ^Jvl-^ l\t(x,it) (xo{t) + f 



2 '^l,! 



n/3 



dt2 Xg(ti)Xg(t2) Xo(il) + 777 2:0(^2) + — + 



n/3 



(3.18) 

This expansion can be thought of as an expansion around slowly varying Ao- One obtains 
a similar term for the anti-fundamental Wilson loop and together they give a one-loop 
effective action which is written in terms of bilinears of V{„.} ( |3.16| ). 
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It is not clear at first sight how to rewrite the comphcated functions V^^.j ( |3.15|) in 
terms of the Wilson loops Wn{x) ( |2.1C1| ), which is a necessary step for constructing the 
expansion (|1.3|), although it is clear that up to global issues the W^s exhaust the degrees 
of freedom of the 9a out of which the ^{n^} are constructed. It turns out that one can 
construct an expansion of V^^.} in a power series in Wn- To leading order in Wn the results 
are simple; in fact one can show that for /c > 2 in ( |3.15| ), 

Vji:l..,^, = 0{W% (3.19) 

whereas for k=l we have: 



1 / d 



ni 



^""'=2™Uj ^'^(^)^OiW'). (3.20) 

The terms of order I > 2 above are of course calculable as well. We will not derive 
( p.l9|) , (|3.20|) here. To illustrate the fiavor of the arguments, we consider the simplest 



non-trivial case of V^m} with Yli=i = 2. There are in this case only two functions 



a=l 
N 



(3.21) 



a=l 



Clearly (gTOD 



= Z^- {l-KufL^. (3.22) 



To illustrate ( ^^201) we have to show that L^ = 0{yV'^). 
To establish that, consider 



,6=1 k-- 



-oo 



2(L„ + L,W^) + 2 [LkWr.->^ + ^^2kl^_j^^ W',WU]- 



(3.23) 



All the Ln are expected to be of the same order in W; we see from (|3.23|) that this order 
must be 0{W'^) and that: 
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Similar considerations allow one to prove (|3.19|) , (|3.20 ) for all 



Hi 



The importance of the observations ( ^.191 ), ( |3.2UD is that to select terms in 



which are bilinear in W (and thus contribute to G*^^-* ( |1.3[ )), we need to keep only terms 
with k = 1 (|3.20|) i. In terms of the expansion in 9 (see ( |3.18|) ) this implies that one need 
only take terms up to first order in 9 from each fundamental trace (or up to second order 
in 9 if the trace is done in the adjoint representation). Higher order (in Wn) contributions 
to the effective action (|1.3| ) can be derived systematically by using ( |3.23| ), ( |3.24| ) and 
generalizations to expand the V|„.} in a power series in Wn- The fact that the winding 
n contribution to the one-loop effective action (|3.10|) is directly related to Gn^ in ( [T7^ ) 



(2) 

together with the above observations allow one to calculate Gn to this order. We now 
turn to this calculation. 

(2) 

3.2. Calculation of the inverse propagator for Wilson loops Gn 

Starting from ( |3.17| ), we have to expand in 6*, average over the fluctuating Xq and 
rewrite the results in terms of bilinears in Wn- Using the results of the previous subsection 
it is easy to see that: 

1 1 ^ poo jrp 2 2 

crl-loopr^i ^ _£_!_ \^ \^ / m^T+2nin9ab(x) 

J n= — oo a,b=l 

+ ^^1^ t^t2 (^^~a[x(tl)]^~b[x(t2)] (^Xo(tl) + ^^ (i;o(t2) + ^^^ ^ 

(3.25) 

where we defined 9ab = 9a — 9b and neglected 0{W^) terms. 



Now we need to expand 9[x{ti)] ( p.l4| ) and average over the Xq. Only those terms 
with equal numbers of derivatives in the expansion of the 9a9b contribute after averaging 
over quantum fluctuations: 



C.l-loopr^T 1 1 dT nl£l-m^T+2nine^i,(x) 

) n= — ooa,b=l ^ 

dt^ dt2 ij2 ^—^0(J\x)9i;\x)Git^ - t2)\G{tl - t2) + 



(3.26) 

where both derivative in G are with respects to ti. Note that the n = contribution to 
the sum vanishes, since ^6'a(a;) = (mod 1). 



Except for the n = contribution which will be discussed separately later. 
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We can shift the variables of integration to make the integral independent of t2 due to 
the periodicity of the propagators G. Doing the remaining integral over ti gives a factor 
of T 



; nil . 
{2i+iy.- 



N 



11 ^ ^ dT r.H^ 



1 + 



/32 



(3.27) 



Next we change variables to w = where a = \n\(3m. The last equation can easily 

be written in terms of the Wn by using (|3.20D . For convenience we write the answer in 
momentum space (an integral over momenta is implied): 

5-'°op[^] ^-I^Y: C-^) ' r ^e-*(-+i)W^.(fc)^_.(-fc) 



(4-. ,.^u ^228) 
{-lyn fPauV u ' 
(2Z + 1)! V2m2 y' ^ 

V. I — 1 

The integral over u can be expressed in terms of modified Bessel functions using the 
defining equation 

1 

2 Jo 

From this we obtain 



^ °° (' 1 V/l / \ ^ -1 



(3.29) 



As mentioned above, to complete the calculation we have to evaluate the zero winding, 
n = 0, contribution. For constant 6a the zero winding sector gives rise to an unimportant 
temperature independent constant, which is usually dropped. In general, however, the 
n = contribution is non - trivial. Indeed, to quadratic order in Wn one finds: 

(3.30) 
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where half of the contribution is due to the fundamental and half to the anti-fundamental 
representations. 

We need to write this in terms of non-zero winding Wilson loops Wn, which can be 
achieved by utilizing the derivation of ( p.24|) and generalizations. In momentum space, 
one finds: 

^ S ^W-„(*)»^_„(-*) f ^e-'- ± (3.31) 

where we have dropped a temperature independent constant, and as usual a momentum 
integral is implied. 

Finally, doing the T integral gives 



(3.32) 

We can now combine (|3.29| ) with (|3.32| ) and find the one-loop contribution to the 



quadratic action for Wilson loops: 



I 



^ i-iyn / k^\ \fa\i-^-\ ,,,,,, , 2 ^„ d j\ 



(3.33) 

(2) 

Adding this to the term in the effective action coming from the classical action we get Gn 
to one-loop order: 

d d 

1 k'^ ( m? \ ^ / / 2 ^ 2 



a J 2 

(3.34) 

Recall that a = m(3\n\. Eqns. ( |3.33| ), ( |3.34| ) are the main result of this section. We now 
turn to the study of some of their properties. 
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3.3. Properties of G^n\k) 

Despite appearances, the limit m of ( |3.34|) is actually regular. To study this limit 
it is convenient to rewrite the sum over Bessel functions in terms of a compact integral by 



reversing the order of the sum and integral in (|3.27|) . This gives 



m 



(47r)t a' 2g'N 



m 



+ 



du 



u 



1+^ 



where for x > 



and for x < 



fix) 



a'^k'^u 



Kd{a) 



2m? 



- r*(e-4'"'»-i'-e-4"=»)/(^^^)) 



X 



1 M 

e 4 



dt — 64 







|x| 



dt—e~^. 
2 



(3.35) 
(3.36) 

(3.37) 



To study the convergence properties of the integrals in (|3.35| ) we have to examine the 
asymptotic behavior of f{x). First note that in both cases f{x) ~ — x as x — ^> 0. Thus, for 
d < 4 the integrals in (|3.35|) are convergent in the UV (for small u). The second integral 
exhibits a standard logarithmic UV divergence which corresponds to coupling constant 
renormalization in d = 4. For x > 0, it is easy to see that /(x) goes to a constant (—2) as 
X ^ oo. Hence for k^ > the integral representation for G*^^^ ( |3.35| ) is convergent even for 
m=0 (and fixed /3) , for which the exponential IR suppression of the integrals (provided by 
the mass) is absent. This is to be contrasted with the infrared divergent results one finds 
if one expands / in a power series in ^k'^u (derivative expansion) and integrates term by 
term. 

oo the integral in ( |3.37| ) converges to 2^/^x. This implies 
> — oo. Putting this asymptotic form for /(x) into ( |3.35| ) 



In the case x < as Ixl 



e 4 



as X 



that /(x) grows as 
for /c^ < and large u gives 



a^\k\ 



du 
^TFexp 

U 2 



1 - 



4m^ / 



(3.38) 



For /c^ < — 4m^ this integral diverges. By rescaling the integration variable we see that in 
this range G^"^^ develops a branch cut: G*^^) oc {4m? — |/cp)('^+^)/^. We will discuss this cut 
below and will see that it has some important consequences. 
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As mentioned above, we find a regular expression for 

as m — > 0. Tliis is easiest 
to see from ( |3.35|) , in the limit a — 0, ajm = (3\n\ fixed. We find: 



1 A;2 / 1 \ 2 



Jo W "'"2 Jo U2 / 



By using the asymptotic behavior of (|3.36|) , ( p.37|) we see that G'n\k^ > 0) is regular. 



The branch cut mentioned above starts now at /c^ = 0. 

Before going on to apply these results to physical problems, in the next subsections 
we briefiy comment on similar calculations for gauge fields and fermionic matter. 

3.4- Fermionic matter 

For fermionic matter we start from (|3.5|) and remembering that the fermions have 
anti-periodic boundary conditions, i.e. po = ^{2n + 1), we derive the analog of ( |3.10|) . 
Since the 7 matrices are anti-commuting operators, it is natural in this first quantized 
approach to introduce world line fermions to represent them. The easiest way to do this 
is to implement a super symmetric generalization of (|3.8|) , introducing superpartners ip^{t) 
for the x^(t). For more detail see [ pT| . 

For adjoint fermionic matter in two dimensions we obtain 



1 /"°° ftT r 



exp 



dt {^-x^ + + 2zi;'Foi^' + ^ - iAoMt)] {^oit) + ^) ) " 

(3.40) 

where Foi = diAo- 

We need to expand to second order in Aq and then average over the quantum fluctu- 
ations of the path, as well as the i(; fields. Note that the i(; fields satisfy 

{VitlWit2))=-GF{ti-t2)9>''' 

(3.41) 

-Grit + T, t') = Grit, t') =sgn(t - t'). 
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After some algebra we obtain 

4 (47r)2 Jo T'+2 

(^1 - £ dh dt2 ^-^iA^'Hx)fmi - t2)t\-lGF{h - t2f 

+ G(ti-t2)(G(tl-t2) + ^)))). 

(3.42) 

The rest of the analysis follows straightforwardly from the bosonic case. The term 
introduced by the world sheet fermions does not effect the position of the cut in the effective 
action. 

3.5. Gauge fields 

We write here only the one-loop effective action obtained from integrating out gauge 



fields in four dimensions. The interested reader may consult the work of Strassler |22] for 
a thorough discussion on how to handle vector fields in this approach. The effective action 
is given by 



11 ^ ^ r dT 

-e 



n= — OO 

(i - dt, J^' dt2 (f; ^-^{A^'\x)r[G{t^ - t2)r' (3.43) 

{-4l + G{h-t2){G{h-t2) + ^)) 

The rest of the analysis is similar to the scalar case with the cut starting aX — as 
expected. 



4. Consequences for the high temperature continuation of the confining phase 

In reference 0, J. Polchinski has proposed to use the form of G^^) to study the 
confining phase of QCD. The idea, inspired by string theory, is to view the perturbative 
calculation of G^"^^ , G*-^-* , . . . which is in principle only valid at high temperature, in the 
deconfined phase, as an analytic continuation of the confined phase to high temperatures. 
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In string theory such analytic continuations are routine; the analogs of the Wilson lines 
Wn, which are winding modes around compact time are governed by an action: 



1 

9st 



Cstr = — M^n(k)M^_n(-k)[k2 + M^{f3)] + 0{W^y, MliP) = (3\^ - C (4.1) 



with C a positive constant. There is no difficulty in formally extending this formula to 

It has been pointed out that calculations such as those of 0, in which G^"^^ is 
calculated only at zero momentum may suffer from IR divergences, since to find S'efi in that 
case, massless particles are sometimes integrated out, and the effective action is in principle 
non-local already at ^ 0. Of course, in four dimensional gauge theory the spatial gauge 
fields that are integrated out are not really massless but rather develop a magnetic mass 
m ^ g'^{(3)/(3 (which is unfortunately incalculable in perturbation theory), and in adjoint 
QCD2 one may turn on by hand a mass for the matter fields (whose role is precisely to 
mimic the above magnetic mass). One would however like to understand whether the 
results of 0, [0 are reliable in the massless limit, and more importantly whether one can 
indeed learn about the confining phase from such perturbative calculations. 

The results of section 3 help to resolve both issues. Regarding the IR divergences, we 
see that they are harmless, even for massless adjoint matter (gluons) in 2d (4(i). Consider 
for example the case of massless, bosonic adjoint QCD2- The inverse propagator 
vanishes when ( |3.39| ): 



1 ^= 



+ 1 r ^e-^ f{e(3'n\) - r — (e-^ - l)/(fcW«) = (4.2) 
2 Jo Jo u 



(see section 3 for definition of / and derivation). In [Q, the last two terms on the r.h.s. 
of (|4.2|) were neglected, because they are formally small at high temperature (of higher 
order in P^g"^) when the sum of the first two vanishes, i.e. when = Ag'^N . We saw in 
section 3 that the IR singularity due to integrating out a massless field results in a cut to 
the left of k"^ = 0. This cut is not dangerous at Ag'^N = /c^ > 0. The full inverse propagator 
( [1.2|) vanishes at k"^ = 4:g'^N[l -\-0{g'^l3'^y\ (up to logarithmic corrections), approaching the 
result of 0, as the temperature goes to infinity. To put it differently, integrating out 
the massless constituents does not introduce infrared subtleties since the Wilson loops Wn 
are tachyonic at high temperature. Thus the calculations of and are valid, contrary 
to the recent claims . 
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Unfortunately, our result ( 3.33 ) seems, at least naively, to invalidate the basic idea of 
using perturbation theory for G^''^ to study the confining phase of QCD. Our calculations 
give a (one loop) glimpse of the analytic structure of the effective action (|1.3| ) in momentum 
space. Indeed, one of the most important qualitative differences between the confining and 
deconfined phases of QCD is the analyticity properties of the Green's functions, such as 
G^'^\k). In the confining, low temperature phase, G^^-* is expected to be an analytic 
function of k"^ to leading order in (compare to ( [4.1|) ), since any singularities would 
have to be interpreted in terms of interactions of the among themselves and with other 
singlet bound states and would be down by powers of the string coupling In the 

high temperature phase, one expects the structure to change drastically. The Wn are no 
longer the natural degrees of freedom, and in particular there are non-singlet operators 
(quarks, gluons) that may couple to W^. Thus one expects G'^^-' to contain branch cuts 
corresponding to pair production of such non-singlet degrees of freedom and to have a 
complicated analytic structure typical of the deconfined phase already in leading order in 
1/N. 

To determine whether the perturbative calculations of 0, correspond to an ana- 
lytic continuation from the confining phase or to properties of the deconfined phase one has 
to study the analytic structure of G^^^ and in particular look for branch cut singularities 
signaling the propagation of constituents. The cuts we found in section 3 are precisely of 
this kind; they seem to correspond to coupling of Wn to two constituents (quarks, glu- 
ons). Thus, the natural conclusion from our analysis is that the perturbative calculations 
performed in 0, and section 3 give the effective action for Aq in the deconfined phase, 
written in terms of the variables Wn, and not, as one would hope, the analytic continuation 
of the confining phasell. The latter would have a very different analytic structure than the 
former, and than what we find (section 3). One can not easily infer any properties of the 
confining phase (other than it being unstable at high temperature) from these calculations. 
In particular the number of degrees of freedom of the string does not seem to be easily 
extractable. There does not seem to be a simple way to calculate in the confining phase 
without actually following G^"^^ up to the phase transition at where all the non-singlet 
singularities should disappear. This requires a non - perturbative analysis, which may be 
feasible in certain toy models fl^ . 



^ In principle one has to investigate the analytic structure of higher order (in g) corrections 
to ( |1.3| ), but it is clear physically that the qualitative picture presented here should persist to all 
orders. 
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5. Consequences for the structure of domain walls 

In the deconfined phase it is natural to express the action in terms of the 9 variables 
( P^) as opposed to Wilson loops. For slowly varying 9a{x), the Lagrangian ( |1.3|) takes the 
form (to one loop order): 



o N 



^ = UE(V^«)' + ^efr(^) (5.1) 

^ a=l 



where the kinetic term comes from the classical Lagrangian ( |3.6| ) and the one loop potential 
Ves is ( p.ll| ), which in four dimensions is (for < 9ab < 1): 



N oo 



(5.2) 

27r^ 

B4:{9ah) 



3/34 

a,b=l 

where 9ab = 9a — 9^ and B4^{x) is the Bernoulli polynomial 

B^{x) = - 2x^ + - 1/30 



The Lagrangian ( |5.1| ) is invariant under the Zn symmetry ( p.8|) ^ + k/N . 
This symmetry implies that there are in fact different minima of the potential 14ff, 
corresponding to all 9a = k/N for k = 0, 1, . . . , — 1. In terms of Wilson loops, these 

1 TT7" 27r ITlfc 

mmima correspond to Wn = e n . 

By analogy with spin systems @], it is natural to study domain walls separating 
regions in space corresponding to different vacua of VefM. In particular, one can attempt 
to calculate the interface energy a, defined as the free energy per unit area of the wall. 

At high temperature, when g{(3) is small, one may hope to use semi classical techniques 
to study these walls |1TT[] . To calculate the free energy of a domain wall between regions in 
space corresponding to Wn = 1 and = exp(27rm/A^), say, we have to find a solution of 
the effective action ( pT|) which behaves as 9a{z — oo) = 0, 9a{z oo) = 1/A^ (for all a). 



It can be shown [|Tl| that the minimal action solution is obtained by choosing a particular 
path in 9 space, parametrized by@: 

9a=qiz)/N, a = l,...,iV-l, 

N-1 (5.3) 

= J^^i^)- 



4 There is some debate in tlie literature regarding the existence and physical significance of 
such wahs p^]- [|l2|, [Q- |16[ |. We will assume that they exist and study some of their features. 
^ Note that in this section is not assumed to be large. 
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In this parameterization the Wilson loop takes the form 

W^n = e^(l + -(e2™^-l)). 

Thus g = corresponds to Wn = 1 while q = 1 corresponds to Wn = e^^. The action 
for q{z) ( ^.31) is (L^ is the transverse area of the domain wall): 

..^^^/..(^(,T.>^a-i*^). (5.4, 

[q] = q mod 1. The solution with the right boundary conditions is 

q{z) = (5.5) 

l + exp(^f^7^//3) 

Plugging it back into (^]^) we find that the interface tension is 

^eff 4(iV-l)7r^ 1 



It is clear from (|5.5|) that the scale of variation of the solution q{z) is ~ 1/g. Hence, 
each derivative comes with a power of g. Formally, this implies that higher derivative terms 
in the effective action (as well as higher loop contributions) modify the solution ( [5 .51 ) and 
a ( p.6| ) only slightly, as the effective coupling at high temperature g{P) is small. One might 
worry [|16[ that infrared effects due to integrating out massless gluons may spoil the formal 



power counting. To examine this issue one has to look at higher derivative terms in the 
effective action for q{z). 

Our results are in general insufficient for this task, since we have not calculated 5'eff 
for arbitrary q{z). However, to study the above infrared issues it is enough to consider the 
behavior of the domain wall profile q{z) as z ^ — oo (say), since then q (^) is small, 
q{z) ~ exp(-y/^^), and one can use our results from section 3 to write the effective action 
to order q'^ (but to all orders in the derivative expansion (|3.14|) )i: 

>5eff(?) = a — / - ^ + 



n'p'dl ^(e-^ - l){f{-dlP'n'u) + 1))) g(.) + 0(g^ ■ ■ ■)• 

(5.7) 



The modifications to eq. (5.7) as compared to eq. ( p. 39 ) are due to the difference between the 



scalar ( |3.4|) and gauge (|3.7|) determinants (see section 3), and should not matter for the qualitative 
remarks that follow. 
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We see that the situation is different from that of section 4. Plugging in the asymptotic 
form of (|5.5| ) in ( pJl) we find that /c^ = -dl = -{N/3)g'^/(3'^ < 0, so the terms on the 



second and third hues of (|5.7| ) are not small corrections to the "leading behavior" obtained 
from the first two terms on the r.h.s. Unlike the Wn in section 4, the mass squared of 
q is positive (and equal to the square of the electric mass). Therefore, the branch cut 



discussed in section 3 significantly alters the behavior coming from (|5.1| ). It seems at first 
sight that ( |5.7|) completely eliminates the possibility of the existence of a domain wall, 
since q{z) = exp(a2;) is not a solution for any real a. However, the correct interpretation 
is different. 

One important effect due to higher order contributions in g that we have neglected so 
far is the generation of a "magnetic mass" for the spatial components of the gauge field 
(static magnetic fields are screened). This magnetic mass, which is of order mmag — g'^/P 
is not perturbatively calculable |2^], but its presence alters the picture following from 
( ^.7|) . As we saw in section 3, it shifts the branch cut from = to = —^ttl^.^^ ^ 
— 4(7^//?^. The scale of the solution ( p.5|) of [|Tl|] is on the other hand the "electric mass" 
rriei — g/ (3. Since g is small in high temperature QCD, mmag ^ "^ei, so that even with the 
modification of (|5.7|) due to the magnetic mass, the solution ( p.5|) is not infrared stable. 
However, as is clear from (|3.35|) (with m = m^ag, and an appropriate generalization for the 
gauge case), there is now a solution q{z) which behaves asymptotically as q{z) ~ exp{az) 
with a ^ rrimag- 

The main question now is whether the scale of variation of the domain wall solution 



remains mmag for moderate q as well. In [|lq] it has been argued that the answer to this 
question is negative since the infrared scale is determined by the mass of the space - like 



gluons Aab which for generic 6a is y ?T^mag + (^)^- that case, the form of the domain 
wall solution q{z) ( |5.5D is only significantly modified at the tails \z\ oo, and the leading 
behavior of the interface tension ( |5.6| ) is not effected [|TB[. However, we believe that the 
scale of variation of the solution q{z) is of order mmag throughout the wall. The point is 
that there are many physical states whose masses are of order m^ag for generic 9 (or q). 
Examples include the space - like gluons Aab with a = b, and gauge invariant combinations 
like TrF^. These can be pair produced at higher orders in the loop expansion and lead, 
as explained above, to a domain wall profile which varies on the scale m ~ mmag for all q. 

Thus, we conclude that infrared effects change the scale of the domain wall solution of 
TlJI from mei (pT|), ( |5.5[ ) to m ~ mmag- Accordingly, the interface tension a, (|5.6| ) changes 



from a ~ to a ^ This can be easiest seen by computing the contribution of the 
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potential term to the effective action ( p.4|) , using the fact that q ~ q{mz) and m ^ ?7imag- 
Higher terms in the effective action give subleading contributions to the interface tension. 
The free energy of the domain wall seems to be much larger than previously believed, and 
is in accord with the expected non - perturbative behavior of asymptotically free gauge 
theory. 

To actually prove the above assertion one would need to derive the effective action 
for finite q and verify that it admits a finite action domain wall solution with the above 
described asymptotics. Due to the non - perturbative nature of m^a.g and other problems 
this seems difficult at present. 



6. Conclusion 

The main purpose of this paper was to study the properties of the effective action for 
Wilson loops winding around compact time in different finite temperature gauge theories. 
The main results obtained are: 

1) We have calculated the quadratic term in the effective action for the Wilson loop 
Wn to one loop order in the gauge coupling constant, and have outlined the calculation 
of higher order (in W) terms. We found that the inverse propagator of Wilson loops G*-^-* 
( p.34|) contained a branch cut which was interpreted as due to pair production of 



constituents in the external Wn field. The improved understanding of the dynamics of 
Wilson loops was then used to reconsider two recent proposals to apply high temperature 
perturbation theory to different physical problems. 

2) We have discussed the idea ^ (see also 0) that one can use perturbative tech- 
niques in QCD to deduce properties of the confining phase, "analytically continued" to 
high temperature. We argued that the analytic structure of the quadratic term in the 
effective action for Wilson loops that we found does not support such an interpretation 
of the perturbative calculations. The inverse propagator in the confining phase is not 
expected to exhibit any singularities in momentum space, to leading order in while 
such singularities do appear in the perturbative results. 

3) We found that certain domain walls between different vacua with broken sym- 
metry that were extensively studied in recent literature [|TD[ - , |]T^ - |T^] are modified 



significantly due to infrared effects. In particular, we argued that the free energy of such 
domain walls behaves like l/g"^ as opposed to 1/g as previously believed. 

22 



There are many natural extensions of this work. One would like to extend the results 
obtained here to higher orders in the gauge coupling. In order to study the large decon- 
finement (Hagedorn) transition one needs, as we saw, to calculate the Green's functions 
in ( |1.3| ) to all orders in g. This may be feasible in toy models of lower dimensional Yang 
Mills theory fl^, where one may use extensions of our techniques to calculate the Hage- 
dorn temperature and perhaps to explicitly see the change in the analytic structure of the 
action ( p. .31 ) between the deconfined and confining phases. The physics of the branch cuts 
in the Wilson loop propagator found above, and in particular their relation to ones that 
appear in the true high temperature vacuum (at different in general), also needs to be 
understood much better. 

It would be interesting to understand what are the implications of the larger free 
energy of domain walls found here for the cosmological scenarios described in the literature 
||T0[] , IT^ . Also, a better understanding of the dynamics of Wilson loops may suggest 
ways to study the disintegration of fundamental strings into their purported constituents 
suggested in |^ . In particular, it would be interesting to obtain a similar picture to that 
found here for unified strings above the Hagedorn temperature. 
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